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Abstract 

We uncover a contribution to the NLO/NLL threshold resummed total cross section for top 
quark pair production at hadron colliders, which has not been taken into account in earlier liter- 
ature. We derive this contribution - the difference between the singlet and octet hard (matching) 
coefficients - in exact analytic form. The numerical impact of our findings on the Sudakov re- 
summed cross section turns out to be large, and comparable in size to the current estimates for 
the theoretical uncertainty of the total cross section. A rough estimate points toward a few percent 
decrease of the latter at the LHC. 



1 Introduction 



Improving the theoretical accuracy of the total cross section for top quark pair production is a ma- 
jor goal, given the importance of top physics, the excellent data-taking ability of the Tevatron and 
the imminent start of the LHC. Yet, relatively few theoretical calculations for this observable have 
been done so far. The most important input are the next to leading (NLO) numerical calculations 
of [1-3], where the NLO correction with accuracy better than 1% was derived. The first exact 
calculation of this observable [4] confirmed these results and their estimates of the numerical un- 
certainties. By demonstrating the appearance of a priori unexpected analytic structures, that work 
also clarified why theoretical progress in top production was hampered for so long. While work 
towards the derivation of the NNLO corrections to the total top quark pair production cross sec- 
tion is underway [5-11] more theoretical effort will be needed before the NNLO result becomes 
available for phenomenological analysis. 

In view of the lack of improved fixed order calculations in the last ten years or so, the only 
other source of refinement in the theoretical predictions for the total top quark pair production 
cross section was based on the so-called soft-gluon resummation [12-17] (related phenomeno- 
logical analyses can be found in [18-20]). The basic idea is that our ability to predict certain 
logarithmic terms to all orders in the strong coupling expansion provides insight into the behav- 
ior of the cross section at higher orders. That, in turn, implies better control over the associated 
theoretical uncertainties. 

Heavy flavor pair production at hadron colliders plays a prominent role in the global soft-gluon 
resummation program. The reason is that at least four hard partons are involved in the underlying 
scattering process. As it is well known (see, for example, Ref. [21]) in the case of scattering of four 
or more colored partons the color algebra cannot be made trivial. This, in turn, spoils the simple 
exponentiation picture familiar from processes like Drell-Yan and Deep Inelastic Scattering. 

As was first established in Ref. [22], the one-loop soft anomalous dimension matrix that con- 
trols the non-trivial soft-gluon correlations in this process diagonalizes in the singlet/octet basis in 
the kinematical configuration in question. This is a very important result that is the basis for the 
simplified exponentiation of the next-to-leading soft logarithms (NLL) in this process. 

Utilizing the above results, the exponentiation formula for the total inclusive top quark pair 
production at hadron colliders beyond the leading logs was given in Ref. [17]. The singlet/octet 
diagonalization mentioned above implies that the soft function A is just a sum A = Ai + As of 
two standard Sudakov-type exponents, which separately describe the exponentiation of the sin- 
glet (resp. octet) color channels. Each one of the two exponents is controlled by its own set of 
anomalous dimensions 

rl JV-1 _ 1 / Mm 2 (l-z) 2 Jo 2 \ 

lnA^AO = / dz— / -|-A f y(a s (^))+A 7; i(a s (4m 2 (l -z) 2 ) . (1) 
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The indices ij refer to the partons in the initial state: ij = (qq,gg), I = 1,8 and N is the Mellin 
moment dual to the kinematical variable p = 4m 2 /s (with s being the partonic invariant mass). The 
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Mellin transform is defined as f(N) = f p N ~ f(p)dp. The anomalous dimensions A,-/ = A,- +Aj 
describe soft-collinear initial state radiation. They have a standard expansion in powers of the 
running strong coupling, see e.g. Ref. [17], and are known in QCD through three loops [23,24]. 
Contrary to A ;J , the anomalous dimensions D,yi control wide angle soft radiation and depend both 
on the initial and final states. They are a priori unknown, but one linear combination of Djj i and 
Djj g can be fixed from existing results for the threshold expansion of the total cross-section [1,2,4] 
at NLO. To fix uniquely both anomalous dimensions at the same order of perturbative expansion, 
the authors of Ref. [17] used heuristic arguments, namely, that the Sudakov exponent for the color 
singlet channel is identical to the one known from processes like Drell-Yan and Higgs production. 
With the help of an explicit calculation, in this paper we are able to directly confirm that assumption 
through NLO/NLL. 

Besides the soft Sudakov exponents A,jj, color dependence is also present in the so-called 
hard coefficients that we discuss next. In the singlet/octet basis one expresses the Sudakov total 
cross-section for top quark pair production at hadron colliders within NLL accuracy as 

o] J OT (N)=o ljA (N)+G IJ .*(N), (2) 
where the two terms are given by 

o 1J j(N)=of°F(N)o? jJ A ihl (N). (3) 

The hard coefficients ajj j are process dependent and once the leading order correction o^°™(Af) 
has been factored out, they contain only iV-independent constant terms. As usual, these coefficients 
are extracted from a fixed order calculation. Their derivation through NLO is the main goal of this 
article. 

In Eq. ([3]) we omitted contributions from Coulomb terms, i.e. terms that in p-space behave as 
~ OCs/P*, where (3 = y/l — p is the small velocity of the quark pair. The Coulomb terms represent 
an effect distinct from the soft gluon logs considered in this article. These corrections have been 
analyzed in Ref. [17, 25, 26] with the conclusion that they have an impact only in the immediate 
vicinity of the threshold. We refer to these references for further details. 

Next we explain the origin of the color index / in the hard functions appearing in Eq. ©. The 
easiest way to see why it should be present is to recall the basic factorization property of gauge 
amplitudes [27,28]. Keeping explicit only information about the color, the factorization relation 
for any ^-particle amplitude M reads 

M I = J-S IJ Hj. (4) 

In the equation above 5/y is the soft function mentioned above, I, J are color indices and H is 
the so-called hard function, which is finite. While the structure of the color diagonal jet function 
J and the soft function S can be made quite transparent based on general process-independent 
arguments [29-32], the form of the process dependent hard function Hj can only be obtained 
from a direct, process specific calculation. The matrix structure in Eq. © naturally translates into 
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differential or fully integrated over the phase space cross sections. An explicit example for that 
procedure can be found in Ref. [33]. 

The color dependence of the hard coefficients of- 1 was not available to any of the previous 
studies of soft-gluon resummation for the total inclusive cross section in hadronic collisions. The 
main goal of this article is to complete this gap in the literature by deriving the exact coefficients 
from a dedicated fixed order calculation, thus allowing a consistent NLO/NLL calculation and soft 
gluon resummation for this observable. It is also a prerequisite for any attempt for going beyond 
the current NLL accuracy level. 

In the original Ref. [17] these coefficients were approximated with the numerically known, 
color averaged coefficient taken from the calculations of Refs. [1,2]. Such an approximation is 
formally correct if one restricts oneself only to the resummation of the NLL soft logs since, as far 
as the towers of logs are concerned, these matching coefficients contribute starting from NNLL. 
On the other side, such an approximate choice is also correct to NLO, since by construction it 
reproduces the fixed order NLO results for the color summed cross-section in Eq.©. Nevertheless, 
one expects that the specific choice does have a numerical impact on the resummed cross-section. 
This is easy to see with the help of the following argument: the total Sudakov cross section is 
a linear combination of two all-order exponents (see Eq.©) with coefficients proportional to the 
hard coefficients o^ r Therefore, a modification of the two coefficients in such a way that their 
color- averaged contribution is kept fixed, results in a change of the weight these two exponents 
(singlet/octet) carry. 

The effect of this modification should be much less pronounced at the Tevatron compared to the 
LHC since there the main production mechanism for top-pair production is through light quark pair 
annihilation. It is known that for the qq production through NLO/NLL only color octet contributes, 
and in this case the hard coefficient has always been known with high accuracy. 

One final comment regarding the hard coefficients of t. As was established in Ref. [4] and 

66) 

also investigated in Ref. [26], the "constant" term in the color averaged gg cross-section extracted 
from [1] differs by around 7% from the exact value. It is reasonable to suspect that this correction 
might be quite sizable. We detail our findings in the next section. 



2 Results 

To derive the hard coefficients ajj j we follow the technique described in Ref. [4] . The only re- 
quired modification consists in the need to insert suitable color projection operators in order to 
separate the contributions where the heavy pair is in singlet/octet state. To this end we define the 
singlet state as 

|l,/,7) = -^8y|l), (5) 
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Figure 1: Color indices in a cut graph for top quark pair production, necessary to define the color 
projection onto singlet/octet states Eq. The dashed line represents the cut through the top line 
with color index i, and the anti-top line with color index j. 



where i and j are the color indices of the quark and anti-quark respectively. The projection onto 
the singlet state |1), can now be performed with the simple replacement 



T a 2 rpb\ rpb 2 _ 1 rpa\ T a 2 T b\ T b 2 



(6) 



where the color indices can be understood with the help of Fig. \T\ The remaining contribution is 
simply attributed to the color octet state. 

It is interesting to note that the color separation is well defined only in the vicinity of the thresh- 
old. Further away, both contributions become separately divergent. This is not really surprising, 
since the octet state will have a tendency to attract radiated gluons and hadronize into a singlet. The 
necessity to combine both singlet and octet contributions in order to obtain finite cross sections is 
only visible starting from O (|3 3 ) , and does not affect our discussion of soft-gluon effects. 

Once the color separation has been accomplished, the result needs to be expanded around 
threshold. The resulting expressions (keeping also the Coulomb terms) read 
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The coupling a s is the renormalized MS coupling evaluated at scale fi 2 and running with «y = n { + 1 
active flavors. We follow the definitions and conventions from Ref. [4]. The relation between the 
coupling running with n t and n { + l flavors can also be found there. 

The results in Eqns. (1718191101 ) are in agreement with the ones extracted in Ref. [26] from 
calculations of quarkonium production at hadron colliders [34,35]. 

For applications to soft gluon resummation the above results are also needed in Mellin space. 
One can easily switch between the two representations of the fixed order threshold expansion and 
the relevant formulas can be found, for example, in Ref. [18]. To further simplify the expressions, 
we have effectively absorbed the Euler constant into the soft function by switching to a modified 
Mellin moment N = /Vexp(yE). After performing the Mellin transformation we can extract the 
exact expressions for the hard coefficients G^- 1 as defined in Eq.® by keeping only the non- 
log (TV) terms. The results read 
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The coefficients Eq. (I15I16I17I18I) are the main result from the present work. Due to the lack of 
a singlet contribution through that order in perturbation theory (see however the discussion at the 
end of Sec. [3]), 8 coincides with the known expressions for the color averaged cross section in 
the literature [1,2,4]. On the other hand, the coefficients for the gg reaction are new. They disagree 
with the corresponding coefficients presented in Ref. [18]. Such disagreement is not surprising 
given the fact that the color dependence of the hard coefficients has not been taken into account in 
the attempt made in that reference to exponentiate the NNLL soft logs. 

Next we comment on the properties of the results above. The vanishing of the LO singlet 
contribution in the qq is well known, and is due to the fact that the LO reaction is mediated by 
an s-channel gluon. Since the vanishing is not due to kinematics but due to color effects, one 
expects that at higher orders this property may no longer be true. At NLO the virtual corrections 
are again zero due to the sandwiching of the one-loop amplitude with the projected tree-level 
diagram. On the other hand, the square of the one-gluon real emission diagrams is not identically 
zero. Our direct calculation establishes that the color singlet contribution in the qq reaction at NLO 
is suppressed by a factor of p 2 relative to the color averaged tree-level contributions and is thus 
subleading. Such leading behavior is due to the absence of Coulomb singularities and stronger 
suppression from the three-particle phase space. 

A rather striking feature of the exact color coefficients is that their color dependence is "stan- 
dard", i.e. they are simply polynomials in C> , Ca, etc. This is to be contrasted to the color averaged 
coefficients used in the earlier literature where color factors ~ 1 / (N 2 — 2) appear. 

It is very interesting to try to estimate the size of the numerical effect of the new terms in the 
gg reaction derived here and in Ref. [4]. To that end we calculate the hard corrections and 
compare them to their counterparts from Ref. [17]. In the calculation we take N c = 3, /j 2 = m 2 = 
m top' a s y 1 / = n i + lj ~ 0.108 and we restore the dependence of Je as explained above. We get 
the following results 

o H(BCMN) = 1 + ^ 14 . 39+0(a 2 ); (19) 

<% (BCMN) |c 3 exact = 1 + ^12.04+0 (a 2 ), (20) 
<^,i = 1 + — 9.16 +0(a 2 ), (21) 
< g ,8 = 1 + ^13.19+0 (a 2 ), (22) 

The function c^ BCMN ^ is defined as in Ref. [17]. The function c^ BCMN ^ \c 3 exact has the same 
functional form as the one in Eq. ([T9l , but the value of the constant C3 has been modified from 
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Ct, = 37.23 (as derived from Ref. [1] and as applied in Ref. [17]) to its exact value C3 = 34.88 (as 
derived in Ref. [4]). 

Dividing Eq. (121I22|) by Eq. (fT9l) and recalling Eq. ©, we see that the effect of the color 
dependence in the hard coefficients is a decrease of the singlet (resp. octet) Sudakov cross section 
by 12% (resp. 3%) compared to the ones in Ref. [17]. 

This is a large effect. Indeed, the shifts we observe in the Sudakov factor are as large in 
size as the present conservative estimate [19] of the theoretical uncertainty of the total top pair 
production cross section and are significantly larger than the total cross section uncertainty estimate 
in Ref. [18]. A detailed phenomenological investigation of the results derived here will require a 
dedicated analysis. 

It is also interesting to demonstrate the impact purely due to the numerical uncertainty in the 
constant C3. To that end we consider the shift with respect to the results in Ref. [17]. Dividing 
Eq. (1201) by Eq. (TT9l) one can easily see that its effect is to decrease the hard function a^ g (and thus 
the whole Sudakov cross section) by 5%. This is also a very significant effect given that its origin 
is pure numerics. 

Another view of the impact of our results can be obtained by looking at the cross section 
expanded to NNLO, similarly to what has been done in [1 8] . Ignoring terms coming from Coulomb 
enhancement, and cutting the logarithmic expansion at log 2 [3 (see discussion in Section [3]), the 
result presented in Eq. (21) of Ref. [18] for the NNLO contribution reads 

4? = <5 S' m (P) (4608 log 4 p+1894.91og 3 p- 3.4811 log 2 p + o(log(3)) , (23) 
where the expansion parameter has been taken to be a$/ (4n) 

a„(P) = og»(P) + + (^) 2 °S } + ■ ( 24 ) 

It turns out that the coefficient of log 2 (3 exhibits an accidental cancellation and is in fact given by 

- 14306.9505 + 384C 3 . (25) 

Inserting the exact value of C3 derived in [4], the coefficient of log 2 P in Eq. (l23l) changes to 

-912.35, (26) 

i.e. the change of only 7% in the value of C3 results in a magnification by a factor of about 260 of 
the coefficient of the quadratic log of o gg ($) at NNLO. 

3 Summary and Implications Beyond NLO/NLL 

In this work we demonstrate that separate color singlet / color octet hard (matching) coefficients 
need to be introduced in the Sudakov total top quark pair production cross section. With the 
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help of a dedicated fixed order calculation we derive these coefficients in analytic form. The 
difference between the hard coefficients for the two color states has not been considered in earlier 
soft-gluon resummation literature. We estimate the effect of these new contributions showing that 
they decrease the Sudakov total cross section by 12% in the singlet and by 3% in the octet channel. 
These shifts are large when compared to the current conservative estimate of the uncertainty on the 
total top quark pair production cross section. 

The effect of these new contributions on the total top production cross section will be somewhat 
reduced due to the subtraction of the O(oc s ) contributions from the Sudakov cross section (see 
Ref. [17] for detailed description of the NLO/NLL matching procedure). However, we have also 
demonstrated that the new corrections modify the terms in the Sudakov cross section at order 
(0Cg ) by a large amount. Therefore, a few percent effect on the total top production cross section 
at the LHC can easily be anticipated. The exact size of the impact of our findings can only be 
obtained from a detailed phenomenological analysis. 

With the results derived here and in Ref. [4] the NLO/NLL program for the total top quark pair 
production cross section at hadron colliders is now completed. The precision requirements for this 
crucial observable for the LHC program are very high and mandate improved theoretical precision. 
In principle, to achieve that one has to go beyond the current NLO/NLL level of accuracy. 

Significant progress has already been made towards the direct fixed order calculation of the top- 
pair production cross section at NNLO [4—11]. In the near future, improvements can be expected 
from the careful analysis of the new results reported here and in Ref. [4]. The natural step beyond 
that is to try to promote the resummation formalism to the NNLL level. We discuss in the following 
how this can be done. Before we proceed, we would like to make a comment concerning Ref. [18] 
where such an attempt has already been made. The comparison with the direct exact calculation 
reported here shows that the one-loop hard coefficients used in that reference are incorrect. Since 
this discrepancy starts at the level of NLO/NLL it will clearly also affect their predictions for the 
NNLL terms. For that reason we have excluded the single log terms from our discussion around 
Eq. d23l). 

As we emphasized in our previous discussion, the basis [22] of the threshold exponentiation 
for the top pair cross section is the singlet/octed diagonalization of the soft anomalous dimension 
matrix in this special kinematics. Therefore, the extension of the resummation formalism to the 
NNLL level, requires to first verify if the corresponding two-loop massive anomalous dimension 
diagonalizes in a similar manner. As of writing of this paper there exists no such indication in the 
literature. In fact, the only known [36] property of the two-loop massive anomalous dimension 
matrix is that it should differ from the corresponding massless one (known through two-loops 
from Ref. [31, 37]) by terms vanishing in the massless limit as powers of the mass. Clearly this 
information is insufficient to determine its behavior near threshold. 

The next open problem in the NNLL resummation program would then be the derivation of the 
two-loop anomalous dimension Aj,8 appearing in Eq. ©. Arguments about its value are given in 
Ref. [18]. In the present work, based on an NLO fixed order calculation, we make no statement 
about it. That would be a subject for future investigation. 
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The contributions from Coulomb singularities through two-loops are known from other pro- 
cesses, and have been summarized in Ref. [18]. 

Finally, we turn our attention to the hard coefficients ajj l from Eq. ©. The complete set 
at NLO has been presented in this work. The corresponding two-loop corrections can only be 
extracted from a future two-loop calculation of the top-production cross-section near threshold. 
Clearly, this is a very demanding task. Moreover, there might be an a priori non-vanishing con- 
tribution from the square of the one-loop virtual diagrams to a q q \ starting from order oc^ (due 
to Coulomb enhancements and weaker suppression from the two-particle phase space). If indeed 
nonzero, it will contribute to a tower of NNLL soft logs and might have a numerical impact on the 
Tevatron predictions. Such a possibility has not been investigated so far in the literature. 

As a final comment we would like to point to our discussion in Ref. [4], where we have argued 
about the rather limited phenomenological value of truncating the all order exponentiation to derive 
partial NNLO (or higher) terms. Specific examples can be found in Fig. 3 of Ref. [38] and Ref. [4]. 
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